Tunneling spin current and spin diode behavior in bilayer system 
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The coherent tunnehng spin current in the bilayer system with spin-orbit coupling is investigated. 
Based on the continuity-like equations, we discuss the definition of the tunneling current and show 
that the overlaps between wavefunctions for different layers contribute to the tunneling current. We 
study the linear response of the tunneling spin current to an in-plane electric field in the presence 
of nonmagnetic impurities. The tunneling spin conductivity we obtained presents a feature asym- 
metrical with respect to the gate voltage when the strengthes of impurity potentials are different in 
each layer. 
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I. INTRODUCTION 

The study of the tunneling process has a long history 
ever since the foundation of quantum mechanics for there 
are many applicable effects. As an example, the coherent 
tunneling of cooper pairs between weakly connected su- 
perconductors, known as the Josephson effect, has been 
employed in the design of superconductor circuits that 
are prospective in quantum computation and quantum 
information [l|, y] . For the future information process- 
ing and storage technologies, the emerging field spintron- 
ics [3|, y, |5| aims mainly at coherent manipulation of spin 
degree of freedom and controllable spin transport in solid 
states. Owning to these anticipations, some attention 
has been absorbed in the study of the coherent spin tun- 
nelingprocess either in the conventional Josephson junc- 
tions 0,0,13] or in the ferromagnetic tunnel junctions [3|. 
However, as we are aware, the junction constituted by a 
bilayer two-dimensional electron gases has not been con- 
sidered yet although there exist versatile features in such 
systems. For example, a resonantly enhanced tunneling 
was reported in the bilayer quantum Hall system when 
the layers are sufhciently close to each other 10]. More- 
over, the spin Hall conductivity, which vanishes for ar- 
bitrarily small concentration of nonmagnetic impurities 
in monolayer, was shown to have a magnification effect 
in the bilayer electron system [ll|]. We therefore inves- 
tigate the coherent tunneling spin current (TSC) in the 
bilayer system with spin-orbit coupling in this paper. In 
our study, the tunneling is included in the unperturbed 
Hamiltonian unlike the conventional tunneling Hamilto- 
nian approach [12| where it is treated as a perturbation. 

The present paper is organized as follows. In Sec.|TTl we 
revisit the definition of the tunneling charge current and 
then give the definition of the TSC with the help of the 
continuity-like equations in bilayer systems. In Sec. IIIIl 
we consider a twin layer system in which the strengthes of 
impurity potentials in each layer are identical. We study 
the linear response of the TSC to the in-plane electric 
field and obtain the tunneling spin conductivity which 
exhibits sharp cusps. In Sec. IIVI the influence on the 
tunneling spin conductivity caused by the variation of 



the strengthes of impurity potentials in different layer 
is considered. The expression of the Green's function is 
extended so as to incorporate such influences. We find 
that the difference between the strengthes of impurity 
potentials in the two layers gives rise to the asymmet- 
rical feature of the tunneling spin current with respect 
to the gate voltage. Finally, a brief summary is given in 
Sec. |V] and some concrete expressions are written out in 
the Appendix. 



II. DEFINITION OF COHERENT TUNNELING 
SPIN CURRENT 

In order to properly define the coherent TSC, let us 
first recall the definition of the tunneling charge current 
across a junction. Such a current is generated by elec- 
trons tunneling from one side of the junction to the other 
side due to the imbalance of the chemical potentials pro- 
duced by an applied gate voltage V. The bilayer system 
considered in this paper resembles the tunnel junction 
and its Hamiltonian is given by 



H = ^I + eVT, + I3t,, 
Zm 



(1) 



where /? is the tunneling strength, / and r^ (a = x, y, z) 
denote the unit matrix and Pauli matrices in the layer 
representation, respectively. The wavefunction of the sys- 
tem is expressed as ^ = (■i/'f , V'b)"^- From the Schrodinger 
equation, we can obtain the continuity-like equations for 
the densities p^ = ''Pe'^ej namely, 



dt dxi 
dt dxi 



— (V-^V/ -h.c), 



h 

if] 



(V'JV'b - h.c), 



(2) 



with j^- = -^^(V'l^V'f - h.c.) and h.c. refers to 
the hermitian conjugation since ^^ should be in a two- 
component form tpg = {ipgpil'ii)'^ if the spin degree of 
freedom is taken into account. Throughout this paper, 



the subscripts i represent the components of a quantity 
in the spatial space and £ label the quantities of the front 
layer with £ — f or the back layer with £ — b. The nonva- 
nishing terms on the right hand sides of Eqs. ^, which 
are caused by the tunneling between layers, indicate the 
overlap between the wavefunctions for different layers. 

Integrating the first equation of Eqs. ^ over the do- 
main of the front layer Vf , we have 

dNf r id f :, 

^T = -^=JjfzdA-j-J ii^li^f. - h.c.) dV, (3) 

where A denotes the area of the layer and Nr = J^ p, dV 

the electron number in the front layer. Here we focus on 
the tunneling-relevant direction, saying the z-direction, 
and assume that there is no in-plane current flowing out 
of each layer. Equation ^ demonstrates that the rate 
of the change of the electron number in the front layer, 
defined as the tunneling current, contains not only the 
conventional contribution Jjf^ dA, but also the overlap 
between the wavefunctions for different layers. The or- 
thogonality and completeness of the states on each side 
of the junction were discussed [ij, [Ij, [l^. Then Eq. ([3]) 
manifests, from another point of view, why the tunnel- 
ing current is conventionally evaluated by the rate of the 
change of the electron number on one side of a junction. 
Now we are in the position to consider the TSC in the 
bilayer system with Rashba-type spin-orbit coupling. A 
natural definition of the TSC density in a system should 
be the rate of change of the spin density Si — ipl s "ipi in 
the £-layer where s = |ct are the spin operators. Here and 
hereafter the overhead arrow represents that the quantity 
is a vector in the spin space. The continuity-like equa- 
tions for the spin density Si in the system with SU(2) 
gauge potentials A i^ an d Ai [l6[ have been obtained in 
our previous paper [ll| . which read 

[jT, - VAfax)Sf + (77- + VAf^x)Jf^ 



'dt 



"dxi 



-T-{iPb s V'/ -h.c), 



(^ - vAwx)Sb + (^ + vAmx)!,, 

= 'A[^}si^,-h.c.), (4) 

with J^j — Re(i/'] \{vi^ s\^() where Vi is the velocity op- 
erator and the curl brackets denote the anti-commutation 
relation. Let ai and ^2 stand for the spin-orbit coupling 
strengthes in the front and back layers, respectively, and 
then we have ^/, = 2^(0, ai ,Q), Afy^-^{ai, 0,0), 

A. = 2^(0, a2 ,0), Aby = -|?(a2, ,0) and Af, = 

AfQ = Abz = Abo = with 1] = h. 

Since the tunneling being considered is spin- 
independent, only a gate voltage can not induce a non- 
vanishing TSC. As we known, an in-plane electric field 
is applied to drive the spin Hall current and the basic 



relation between the spin current and the electric field is 
given by [13] 



^i ~ ^aij'^s^jj 



(5) 



where a^ is the spin Hall conductivity and e^j^ the totally 
antisymmetric tensor. It demonstrates that the flow di- 
rection and the spin-polarization direction of the current 
as well as the direction of the electric field are always 
perpendicular to each other. Since the tunneling is with 
respect to the z-dircction and the in-plane electric field 
is along the x-direction, we focus on the component of 
the TSC polarized in the y-direction. From Eq. ^, we 
have 



dt ^dz 



^l^v)Jl - I 



(^s^V'/ 



h.c), (6) 



where we have taken that the nonvanishing components 
of the spin current are J^ = —Jfy as the electric field is 
along the a;-direction. The last term on the right hand 
side of Eq. ([6]) can be regarded as the overlap between 
the eigenfunctions of s^ . Unlike the case of the tunnel- 
ing charge current, the contributions of the states with 
spin parallel and anti-parallel to the y-axis have oppo- 
site signs. The covariant derivative ^ -t- rjA^y is in place 
of the conventional derivative, which indicates that the 
spin precession due to the spin-orbit coupling leads to 
additional contribution to the TSC. 



III. TUNNELING SPIN CURRENT IN 
TWIN-LAYER SYSTEM 

In the previous section, the definition of the TSC in 
bilayer systems with Rashba-type spin-orbit coupling in 
each layer has been discussed by taking account of elec- 
trons' coherent tunne ling . Unlike the tunneling Hamil- 
tonian approach [12, Tsl in the study of the tunneling 
charge current, here we take the tunneling term in the 
unperturbed Hamiltonian. Since the tunneling strength 
is not necessarily weak in our approach, we are able 
to attain more information beyond the perturbation ap- 
proach. In the present approach, the unperturbed Hamil- 
tonian in the second quantization form is given by 



H = 







ail + a_r. 



:)®(V"-fcr.^")}Ck + Mm. (7) 



We have adopted notations a+ = (ai + q;2)/2, a 
(«! — a2)/2 and C^ = 



(.^kJT'''k,/i''^k,bT'''k,6i) ^^^"■ 



being the creation operator for a spin-up (spin- 



el 

k,£T(i) 

down) electron in the ^-layer. Throughout the paper, 
the boldface of a quantity manifests that it is a two- 
dimensional vector, e.g., k = (kx,ky). The last term 
in the Hamiltonian characterizes the interaction between 



electrons and impurities. The influence of nonmag- 
netic [13, [20, UM, [231 and magnetic [23, UM impurities on 
the spin Hall conductivity were discussed in monolayer 
systems. 

For simplicity, we consider the nonmagnetic impurities 
in both layers are all alike in this section. Thus the poten- 
tial energy of impurities is given by Vim = u J2i '^(r ~ R-i) 
with u being the strength and Ri the position of the 
impurity. Here we assume the interaction strength u is 
weak so that the Born approximation [25] is applicable. 

In the chiral representation, the Hamiltonian is diago- 
nalized by a unitary matrix U with eigenenergies 



£1 = 
£3 = 

with e/c 



fife - 

= Cfc 



a+k - A12, 
- a+k - A34, 



£2 = fife - a+k + A12, 
£4 = Cfc + a+k + A34, 



h-'k-' 



Ai 



VW- 



= ^(eT/-a_fc)2+/32 and A34 = 
Hereafter, we set h — I for sim- 



■a-k)^+l3^. 

plicity. The free retarded Green's function in this chiral 
representation is given by GS^^Jkjw) = diag((cj — £1 + 

i77)"\ (w - £2 -I- ir])~'^,(u} - £3 + i?7)"\ (t^ - £4 + «??)"^) 
with irj being an infinitesimal quantity. And the Green's 
function in the original representation is related to that in 
the chiral representation by the unitary transformation, 
namely, G(k,a;) = C/(k) G(,;.)(k,a;) C/t(k). 

As a macroscopic quantity, the TSC is expected not to 
be affected by the details of impurity location {Ri}. We 
assume there is no correlation between impurities and 
employ the impurity averaging techniques and the di- 
agrammatic method |25l . |26[. A physical quantity Q 
for the whole system is obtained by taking average over 
impurities' configuration, namely, Q = (Q({Ri}))im = 
Hi J ^^Q({Ri}). In the Born approximation, the av- 
eraged retarded Green's function G^(k, w) satisfies the 
Dyson equation 



G«(k,a;) = G^(k,c^)-KG^(k,c^)(i 



U^riiyn 



^G«(q,c.))G«(k,c.), (8) 



A 



where nim stands for the impurity concentration. 
The above equation has a self-consistent solution, 
'^(^h)(^'^) = diag(gi,52,53,54) with g^ = (w - e, - 
wnim + ^)^^ where the subscript j runs from 1 to 4. 
Here r = {2Tru'^ni^Np)~^ is the momentum-relaxation 
time and Np the density of states at the Fermi surface. 
Note that the averaged retarded Green's function is diag- 
onal in the chiral representation. We will see in the next 
section that the difference between the strengthes of im- 
purity potentials in two layers results in the emergence 
of off-diagonal elements in G^, > . 

The TSC defined in the previous section can be ex- 
pressed in a symmetric form with respect to both layers, 
namely, the difference between the rate of change of the 



spin operator in each layer 
l.dSl, l,d 



J^-- 



2' dt 



2^dt 



^Gt(T,®s^)Gk). 



(9) 



Its linear response to the external in-plane electric field, 
the tunneling spin conductivity ay , can be obtained by 
using the Kubo formula 

'^{'^) = n 7 / duj iTi {n p{uj + u)i) 

2,1:10 A J 

X [GR{uj + c^i) - G^{u + c^i)]iG^K)j1 



+n^K)[G«K)-G^(c^i)blG^(co + c^i)jj,(10) 

where jf = ^[{a+r^ + a_I) (g) kycr^ — (3Ty (g) ay] and Tr 
refers to the trace taken over the spin indices as well as 
the summation over the momentum. G is the advanced 
Green's function, Up the Fermi distribution function and 
je = evx the charge current operator. 

In the uncrossing approximation [2l| , the dc tunneling 
spin conductivity av at zero temperature is calculated 
as the sum of contributions of the one-loop diagram ctq 
and a series of ladder diagrams a\. The former is given 
by ^ = -2^Tr(G^(k)je(k)G^(k)jl(k)), and denoted 
diagrammatically as 



-T7 1 



27rA 



G-^ 

C3' 



QR 



A direct calculation gives rise to 



tr, 



V _ 
— 



-2-E 



• 2 

sm Lp 



T — ai52 + T — 0354 + a+k{ai -|- a2){gz + 54) 
^12 A34 



X Im 

?2 I /„,2 I „,2 \l2 



0' + {a\ + ai)e - eV^^^^k 
2AI^ 

(3^ + {al + al)k^ + eV^;"' 



(fli -a2)(.g3+54) 



2A34 
{e^V^~l3^-alk^)a+k 



-(fli + a2)(.g3 -54) 



A12A 



12A34 



(ai -02) (53 -94) 



(11) 



where the polar coordinates k — {k cos if ,k sin (f) is 
adopted and ai are the matrix elements of the advanced 
Green's function in the chiral representation. 

In the limit of large Fermi circle, n ^ Aij,l/T with 
Ay = Si — Sj, the vertex correction to the conductivity 
(t| is dominated by the terms with one advanced and one 



retarded Green's function [2l|, which can be expressed 
diagrammatically as 





je + 3 



G^ G^ 

If we introduce a vertex J^^ which is the sum of vertex 
corrections to j|, 



jy 




then a^ can be expressed in a compact form, 
' Tr(G«(k)je(k)G^(k)J'^), 



27rA 



(12) 



and the exphcit expression of cr^ in terms of Green's func- 
tions is given in the Appendix. The momentum inde- 
pendent J^y can be obtained from the transfer matrix 
equation 



J" = 



A 



5]G-4(q)(j^^(q) + :r^)G«(q). (13) 



Figure [3] shows the contour-plot projections of ay 
which are plotted as a function of both the gate volt- 
age and a_. These figures indicate that the resonant 
peaks in a^ begin to converge as a_ decreases. 



-a_ = 5*10" eVm 

_= 1.5*10"" eVm 
-a = 0.45*10"" eVm 




-0.08-0.04 0.04 0.08 -0.08-0.04 0.04 0.08 
gate voltage (mV) 

FIG. 2: (Color online) Tunneling spin conductivities ay are 
plotted as functions of the gate voltage V with different values 
of Q_ given in the legend. The tunneling strength is /? = 
5 X 10~ eV, the momentum relaxation time are r = 660fs in 
panel (a) and r = 66fs in panel (b) while the other parameters 
are the same as those in Fig. [T] 



Here the summation over the momentum can be evalu- 
ated by taking it as integration in the limit of large Fermi 
circle. Then we obtain the tunneling spin conductivity 
ay which is plotted as a function of (3 in Fig. [TJ We can 
see that the TSC vanishes when the tunneling is absent. 



T = 660fe 



= 66fe 





FIG. 3: (Color online) Contour plot of ay which is a function 
of both the gate voltage and q_ . The color bar represents the 

Stt. 



/J(cV) 



FIG. 1: (Color online) ay as a function of /3. The parameters 
are given by F = 0, a+ = 5.5 x 10~^'*eV m, a_ = 0.45 x 
10^^ eV m, ^ = 0.1 eV, r = 660 fs and the effective mass 
m = 0.065 mo as in CaAs with mo being the mass of the free 
electron. 

We plot ay as a function of the applied gate voltage 
V with different a_ and momentum relaxation time in 
Fig. [2 It is shown that a^ is an even function of the gate 
voltage. For small a_, a peak in ay appears around V — 
0. As a- increases, this peak begins to split and its height 
increases. Additionally, the nonmagnetic impurities tend 
to suppress these peaks, as shown in Fig. [DJb) where 
the momentum relaxation time is a tenth part of that in 
Fig. Ha). 



IV. TUNNELING SPIN CURRENT 
ASYMMETRICAL TO GATE VOLTAGE 

In realistic samples, the strengthes of impurity poten- 
tials in those two layers may not happen to be identical. 
We thus introduce Ur and Mj, to denote the strengthes of 
impurity potentials in the front and back layers, respec- 
tively. Accordingly, the interaction between electrons 
and impurities is given by 



Vi„, = ^5(r-R,) 











(14) 



Note that the unit matrix in the layer space is no longer 
appropriate in the present case. This implies that the line 



which refers to the impurity potential in the Feynman di- 
agram (see Fig.|3]) represents a matrix diag(uj/^,M^/v4) 
rather than a number uj A. Accordingly, the Dyson equa- 



A^ Go(k) 




Go(k) 



The transfer matrix equation for the vertex J"^ is also 
modified as 



^ 



^'-^Y.K (t ^ ]^^^'^ 



x(jI(q)+J'^)G«(q)( \ t 



(17) 



FIG. 4; The first order Feynman diagram in the expansion 
of the impurity-averaged Green's function where the momen- 
tum is conserved at the vertex. The vector refers to the free 
Green's function. And the dotted line represents the impu- 
rity potential which is in the matrix form for the unequal 
strengthes of impurity potentials between layers. 



tion for G^ is then written as 
G«(k, uj) = G^ik, Lu) + G^(k, Lu)N 



^ 



E 



A I JG«(q,c.)(^ A I )|G^.(k,^).(15) 



We obtain a self-consistent solution for the above equa- 
tion 






Rii Rf 

Rf R22 

i?33 Rb 

-Rb Rii 



(16) 



with nonvanishing off-diagonal elements even in the chiral 
representation. The explicit expressions for those matrix 
elements are given in the Appendix. The momentum 
relaxation time for the front and back layers are simply 
given by Tj = {2Tru'jni^Np)~'^ and r^ = (27rM^ni,„iVp)"\ 
respectively. 

The difference between the strengthes of impurity po- 
tentials in two layers not only modifies the diagonal el- 
ements of G?^^ but also requires the appearance of off- 
diagonal elements inevitably. Now the contributions to 
ay by the diagonal elements can be directly obtained 
by replacing gi and a^, respectively, by Ra and An in 
Eq. ini) and Eq. (K2^ . Here An = R*- denote the diago- 
nal elements of the averaged advanced Green's function 
in the chiral representation. And the contributions to CTq 
given by the off-diagonal elements reads 



<d = -2/?E^i'^''^I"i 



xi?f (Aii-yl22 ( a+) + {A11+A22) ^— : - 

m A12 

+ (^33 + ^44) ^ — - (^33-^44)^-^) 

-Rh{{A33-A4A){—+a+) - (^33 + ^44) ^^^ 

m A34 

~{Aix+A22) r — + (^11-^22)^^ ) 

A34 A 34 



And the terms in a^j^ given by the off-diagonal elements 
is written out in the Appendix. 




0.04 



^ -0,04 _o,os' 
gate voltage {mV) 



^ a (10-"SVmJ 



o^l-IO^/Sn) 




0.08 0.04 



-0,04 _o,08 
gate vollage (mV) 



a_(iO""eVm) 



FIG. 5: (Color online) ay as function of V and q_. The 
relaxation-time differences are Ar = 8 x lO^'' in panel (a) 
and Ar = 5 x 10"'* in panel (b) with t\, — 660 fs. 

The influence of the difference of impurity potentials 
on the tunneling spin conductivity can be observed in 
Fig. [5l Here we introduce the difference in relaxation 
times At — {tj — Tb)/Tb which is taken to be 8 x 10~^ 
and 5 x 10"** in Fig.[5Ka) and Fig. [SKb), respectively. As 
the difference Ar increases, one of the resonant peak in 
ay tends to be suppressed and finally it becomes a valley. 

Figure [5] shows that the variation of the strengthes of 
impurity potentials between layers leads to the asymmet- 
ric dependence of the TSC on the gate voltage when the 
in-plane driven electric field is fixed. This unilateral con- 
duction feature makes the bilayer system with different 
strengthes of impurity potentials a candidate for the re- 
alization of spin diode. 



V. SUMMARY 

We have studied the coherent TSC in the bilayer sys- 
tem with spin-orbit coupling. We first revisited the def- 
inition of the tunneling charge current with the help of 
continuity- like equations in the bilayer system since the 
tunneling between layers causes the nonconservation of 
the density in each layer. In additional to the conven- 
tional contribution, the tunneling current contains those 



related to the overlap of wavefunctions for different lay- 
ers. This is intuitional for us to define the coherent TSC. 
We showed that the contributions of the wavefunction 
overlaps for states with spin parallel and anti-parallel 
to reference axis have opposite signs. Unlike the con- 
ventional tunneling Hamiltonian approach, the tunnel- 
ing strength in our study is not necessarily week since 
it is treated as a part of the unperturbed Hamiltonian. 
In the light that only a gate voltage can not induce a 
nonvanishing result, we studied TSC in response to an 
in-plane electric field. The spin conductivity was cal- 
culated in terms of Kubo formula by taking account of 
nonmagnetic impurities. We firstly investigated twin- 
layer systems and then considered the effect caused by 
the difference between the strengths of impurity poten- 
tials in different layers. Meanwhile, we developed the 
techniques dealing with the impurity-averaged Green's 
function in the bilayer system. We showed that there 
must exist nonvanishing off-diagonal elements in the av- 
eraged Green's function even in the chiral representation 
if the strengthes of impurity potentials in two layers is 
different. Sharp cusps in the tunneling spin conductivity 



appear near the null voltage and are suppressed by the 
impurities. We found that if the strength of impurity 
potential in one layer is different from that in the other 
layer, the TSC exhibits the asymmetrical feature with re- 
spect to the gate voltage. This reveals that the spin diode 
can also be realized in the bilayer system with different 
strengthes of impurity potentials in different layers. 
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APPENDIX A: EXPRESSIONS FOR SOME 
COEFFICIENTS AND THE MATRICES 

The concrete expressions for the matrix elements of the 
averaged retarded Green's function is given by 



J 



eV—a- k 



Rii = 

R22 = 

R33 = 

i?44 = 

Rf = 

Rh = 









UJ 


-t2 


-w+ 


A12 


—w^ 














(^- 


-ei)(uj- 


-e2)- 


- w+(2ti) — 


ei - 


£2) + 


((jj+ + UJ 
A12 


k 
— w_ 


- LU- 


-)- 


- 2uj- 


-{eV - 


a- 


-kV 


(c.- 


-£i)(c.- 


-e2)- 


- w+(2a; — 

UJ 


£1 — 

-£4 


£2) + 
- W7+ 


A34 


k 
— W^ 


- UJ- 


-)- 


-2uj. 


.{eV - 


a- 


-ky 


(c- 


-e,)icj- 


-e,)- 


- w+(2a; — 


£3 ^ 
-£3 


£4) + 


+ A34 


k 
— W^ 


- UJ- 


-)- 


-2uj- 


-{eV + 


a- 


-ky 


1^- 


-e3){uj - 


-e,)- 


- Wj^{1lo — 


£3- 


£4) + 


(w+ + tJ 


-)i^+- 


- UJ- 


-)- 


- luj- 


-{eV + 


a- 


-ky 












A 


12 

















{(jj — ei){uj — £2) — w+{2uj — £1 — £2) + {uj+ + cj_)(cj+ — uj^) — 2uj^{eV — a^k) ' 



A34 



(uj — e^){to — 64) — w^(2uj — £3 — 64) + (c<j+ + tu-){ujj^ — uj-) — 2uj-{eV -\- a-k) ' 



(Al) 



«th w± — ^{{uf ± Mb)nim — ■:f^) where t± — ^(t/ — Tt) are introduced. 



When the strengthes of impurities are equal in the two layers, the vertex correction to av in terms of Green's 



functions is given by 



T^ 



5;»V{R.[(i- 



a-(eV — a^k)^ ,k a^leV — a^k) 

^- '-)aigi + {- + a+ ^ '- 

Ai2 m Ai2 



)a2g2-a+(ai + a2)(53+34) 



,k a(eV+a^k), ,k a^(eV + a^k), e^V^ + [3^ - a^k^ , 

-(—+"+ ^ jasffa-l — + a+H r^ )a^gA-a+ : — ^ (ai-a2J(53-54J 

m A34 m A34 A12A34 

,eV — a-k , ., . eV + a^k ' 

+a-[ ^ (ai - a2)[g3 + gi) H ^ (ai + a2)(g3 - 04)) 



A12 



A 



34 



i( J12 + J34) 



+Re 



eV — a-kk a-(eV — a-k) k a-(eV — a^k) 

r (( a+ + r jaigi-{ a+ )a252 - a+ (01-02) (.93 +54)) 

A12 m A12 m A12 



eV+a^k k a_{eV+a^k) k a_(eV+Q!_A:) 
I (( \-a+ )a3.g3-( l-a+H ^ )a454 + a+(ai+a2)(53-54)) 

A34 TO A34 TO A34 



n2 ^^34 



«i.92 a3.94^ a-{e^V^ — l3'^-a^k'^) 
Ai 9 A 



A12A, 



12^^34 
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where J12 , J14 , J32 , ^34 are the matrix elements of the vertex J^^ which can be obtained by solving Eq. (I13p . 

The difference of strengthes of impurity potentials in each layer brings the off-diagonal elements to the Green's 
function in the chiral representation. The corresponding contribution to the vertex correction of the conductivity is 
given by 
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